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Introduction. The concept of the sohton and the inverse scattering transform 

approach to solving certain non-linear equations were introduced twenty-five years ago 
in two landmark papers by Zabusky-Kruskal and Gardner- Greene- Kruskal-Miur a [1,2]. 
Since then, soliton theory, or the theory of integrable systems, has had enormous impact 
on applied mathematics and mathematical physics. Water wave theory [3], nonlinear 
optics [4] , field theory [5] , and relativity [6] are but a few of the areas which have been 
influenced by these ideas. An important aspect of the study of soliton equations is that a 
single example can be of interest in a wide variety of contexts; this universality may well 
be related to the fact that such equations frequently have an underlying geometric mean- 
ing. For instance, the sine-Gordon equation, which first appeared in differential geometry 
[7], is also used as a model for dislocation of crystals [8], for field theory [9], and nonlin- 
ear optics (self-induced transparency) [10]. In the same vein, several authors [11], [12], 
[13], [14], [15], [16], have developed the connection between the Serret-Frenet equations 
and other elements of classical differential geometry and various well-known integrable 
models, including the non-linear Schrodinger equation, and the modified Korteweg-de 
Vries equation. 

Continuing this geometric theme, we describe in this paper certain families of 
integrable equations which govern the motion of curves in the plane, on the sphere, and 
in three-dimensional space. We focus our attention on evolution equations which preserve 
distinguished local invariants of the evolving space curve or associated swept-out surface. 
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Our point of departure is the localized induction equation (LIE): •yt = 7s >^ 7ss — 
kB. lie is an idealized local model of the evolution of the centerline of a thin vortex 
tube in a three dimensional inviscid incompressible fluid. Here, the subscripts t and 
s denote partial differentiation with respect to time t and arclength s for the evolving 
curve 7(5, t) in R^, which has curvature k, torsion r, and Frenet frame {T,N,B}; the 
multiplication sign denotes cross product. For derivation and history, we refer to [17], [18]. 
More accurate, non-local models have been considered; we refer the reader to the work of 
Moore-Saffman [19], Klein-Majda [20] and the references mentioned therein. LIE has a 
second interpretation as a "potential" form for the classical one-dimensional continuous 
Heisenberg ferromagnet [21]. 

Hasimoto [22] uncovered the relation of LIE to soliton theory by showing that it 
induces an evolution on the complex curvature ijj = ne^ J ^ " governed by the nonlinear 
Schrodinger equation (NLS): ifjt = i{i^ss + IIV'PV')- NLS is a well known example of 
an equation with soliton solutions. Since Hasimoto's discovery, the structure of LIE has 
been more fully spelled out in the context of inflnite dimensional Hamiltonian systems. 
In fact the Hamiltonian for LIE is just the arclength functional on curves (as shown by 
Marsden-Weinstein [23], who also described the relevant Poisson structure). Further, 
associated to LIE is an infinite sequence of commuting Hamiltonian vector fields, with 
Hamiltonians which can be expressed as global geometric invariants of the curves. All 
of these equations are of the form 7t = VF = aT + bN + cB, where VF is a geometric 
vector field, which means the coefficients a, 6, c are functions of k, r, k' = k^,t' = Ts , and 
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higher derivatives with respect to s. We hst the first few terms of the localized induction 
hierarchy (LIH) [24] of commuting vector fields, as weU as their associated Hamiltonians 
(the vector field Xq is exceptional): 
Xn = -T, 



1 



Xi = kB, /i = / (is, 

.2 



X2 = —T + k'N + ktB, l2= -t ds, 
2 Jj 

X3 = k^tT + {2k't + kt')N + {kt^ - k" - -K^)B, 
-K^ ds, 

7 2 

^ T ' 2 8 ^ 

+ {kt^ - 3(kV)' - -/tV - kt")B, 
I4 = I -K^rds, 



LA 

/ 2 
X^ =7^(X,_l), n>0. 

As the last line suggests, LIH is generated by successively applying a recursion operator 

starting with Xq. For backround material on recursion operators in the theory of inte- 

grable systems, see [25]; the recursion operator for LIH was introduced in [24], and is 

defined below. 

The vector field X2 also arises in a related fluid mechanical context, via a refined 

version of LIE. Fukumoto-Miyazaki [26] derived a "localized induction equation" for thin 

vortex tubes, which allows for axial velocity; their equation be expressed in terms of the 
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vector fields of LIH via 7t = {CX2+X1), c some parameter. Lamb [11] had previously con- 
sidered the (essentially equivalent) evolution equation 74 = X2 — 3roXi + Sto^Xq, which 
has the following special property: starting with an initial curve 7(5, 0) having constant 
torsion tq, the evolution preserves the condition r = constant = Tq, while the curvature 
function K(s,t) satisfies the modified Korteweg-deVries equation Kt = Ksss + ^k^Ks, an- 
other well known soliton equation. Lamb's result can be interpreted as saying that there 
exist distinguished solutions to the Fukumoto-Miyazaki equation which can be completely 
described by tracking just one evolving functional parameter (curvature) rather than the 
usual two (curvature and torsion). The special case tq = corresponding to an evolu- 
tion of planar curves has been discussed in several recent papers [15], [16], [27]; in fact, in 
[28], this planar evolution equation is shown to be a "localized induction approximation" 
associated with vortex patches in ideal two-dimensional fluids. Part of our aim here is 
to demonstrate that Lamb's example is but one of a whole family of evolution equations 
preserving local geometric invariants, that such examples are geometrically significant, 
and that they are intimately related to the recursion operator TZ. 

The Recursion Operator TZ and Variation formulas. Formulas for the 
variations of geometric invariants of a curve 7(5, t) evolving by a vector field 7t = VF have 
been derived by various authors [16], [24], [29], [30]. Much of the present paper depends on 
the remarkably simple expressions for these variations in terms of the recursion operator 
TZ in case W is locally arclength preserving (LAP), i.e., in case < Ws, T >= 0. In addition 
to curvature and torsion, we consider natural curvatures u and v related to n and r via 
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u + iv = ip = Ke J '^^ ". The corresponding natural frame {T, U, V} is related to the 
standard PYenet frame hj U -[-iV = {N + iB)e J '^^^' " and satisfies the natural Frenet 
equations Tg = uU + vV, Ug = —uT, Vg = —vT. We begin our Ust of formulas with the 
definitions of TZ and a parameterization operator V, which takes an arbitrary vector field 
X ^ aT + hN + cB = fT + gU + hV to an LAP vector field: 



(a) V{X) = / {Kh) dsT + bN + cB 

[gu + hv) dsT + gU + hV, 
ih) TZ{W) = -V{TxW'), 

(c) W{k) = <-TZ^{W),N>, 

(d) W{r) = < -n^{W), B/k >', 

(e) W{u) = <-n^{w),u>, 

w{v) =< -n'^{w),v> . 

Note that at each time t, u + iv is only determined up to multiplication by a 
complex unit; the formulas for W{u) and W{v) should therefore read "for some choice of 
natural curvatures u{s,t) and f(s,t)". Also, the appropriate choice of antiderivative in 
the definition of V depends upon the class of curves under consideration. In particular, 
in the asymptotically linear case, it is the "antisymmmetrized" antiderivative operator 
J^ f ds — F{s) — |(F(oo) + F(— oo)), F' = /, which gives the terms of LIH hsted above. 
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The formulas stated above were derived in [24] and [31], where similar formulas were also 
derived for the evolution of the frame {T, U, V}. 

Curve invariants. We now present the integrable hierarchies Xn,Yn, Zn,^n, 
UnjCJrn^m of vector fields preserving geometric invariants of the evolving space curve 
7. As noted below, the X^ and Y^ hierarchies have been discussed previously in the 
literature (we include them for completeness); the other hierarchies are new. 

(i) Locally arclength preserving: For n > 0, the LIE hierarchy X^ is locally ar- 
clength preserving. The X^ satisfy the recursion relation X„_|_i = 7^(X^), starting with 
Xq. li W = Xn, then the induced evolution on complex curvature ifjt = W{'iIj) is the 
corresponding element of the NLS hierarchy; in particular, Xi — kB induces the the 
NLS equation (details appear in [24]). Obviously, any linear combination of the X^ will 
also be locally arclength preserving. Our hierarchies will mainly be of this form. 

(ii) planarity preserving: For n > 0, the even fiows in LIH Yn — X2n are planarity 
preserving (the evolution of a planar curve stays planar). The y— recursion operator 
is just 7l'^{W) = —V{Wss) restricted to planar vector fields along planar curves. The 
starting vector field in the sequence is Iq = ^o- ^^ W = Y^ then the induced evolution 
on curvature Kt = W{k) is the corresponding element of the (mKdV) hierarchy; in 
particular, Yi induces the (mKdV) equation itself [15], [16], [27]. Observe that planarity 
preserving is equivalent to preserving the condition r = 0. 

(Hi) constant torsion preserving: For n > 0, the vector fields 

2n+l 



fc=0 ^ ^ 
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preserve the constant torsion condition r = tq. Each Z^ contains a "phantom" or formal 
vector field component X_i; by definition, TtiX-x) = Xq. X-i has trivial geometric 
action; it is introduced for algebraic purposes, to facilitate the recursion formulas for the 
Zn- The Z-recursion operator is {IZ — tq)^; the starting vector field in the sequence is 
Zq = Xq — tqX^i. We remark that along constant torsion curves, the operator IZ — tq, 
thought of as an integro-differential operator on vector fields W = aT + bN + cB, has 
coefficients independent of r. li W — Z^ then the induced evolution on curvature 
Kt = W{k) is the corresponding element of the (mKdV) hierarchy; in particular Zi 
induces the (mKdV) equation, and we recover the result of Lamb [11]. When tq = 0, we 
are in case (ii). 

A closely related sequence of vector fields is given by the "planar-like" vector fields 

2n-l 



^"= E f^^fc ^y-rofX,, 



k=0 

SO called because they have no binormal component and their coefficients are independent 
of r; in fact, the Qn look identical to the Y^, the only difference being that they are 
defined along space curves rather than planar curves. We have O^n+i = Ti^Zn- Again, 
the recursion operator is (TZ — tq)^ and the starting vector field in the sequence is Oi = 

2 

X2 — tqXi = ^T + k'N. The form of fin, together with formulas (c) and (d), imply 
that the Z^ are constant torsion preserving and induce the (mKdV) hierarchy on the 
curvature k. 

(iv) Constant natural curvature and sphericity preserving: Note that if a natural 
curvature is constant {v — vq along a curve 7) then the natural Frenet equations imply 
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c = 7 + {l/v)V is a constant vector, so 7 lies on a sphere of radius 1/v centered at c. 
(Compare this with the standard sphericity condition k~'^ + {k'Yt~'^k~^ = constant, 
given in most differential geometry books, e.g., [32]). Conversely, if 7 lies on a sphere 
of of radius r, then a natural frame {T, U, V} along 7 is readily obtained by letting U 
be normal to T and tangent to the sphere (and so V = T x U is normal to the sphere); 
then the corresponding natural curvature v is constant, t) = 1/r, and u is the geodesic 
curvature of 7, regarded as a spherical curve. 

We define a constant f -preserving analogue of the planar hierarchy Y„,: cui = 

2 
^T + UgU, etc., are obtained by simply replacing nhy u and N hy U in the expressions 

for the corresponding Yn- From this description, it is clear that the Un define evolutions on 

spherical curves; further, they are generated by an intrinsic recursion operator, S{W) — 

— P(V|i(VF)) - here V is covariant differentiation on the sphere - since Vt(7") = uU , 

and VriU) = ~uT look just like the planar Frenet equations. In terms of TZ, we have 

S{W) = ilZ'^ + v'^){W)] specifically, if we consider the class of asymptotically geodesic 

curves - u decays rapidly to zero - the antiderivative operators in both formulas for S 

are the "antisymmetrized" ones, mentioned above. 

We note that, just as Euclidean elastic curves give solitons for LIE evolving by 

screw motion, "spherical elastic curves" [33] simply rotate when evolving via 7^ = wi. For 

a general spherical curve, one might expect u to evolve according to mKdV, by analogy 

with the planar case. This is not quite true; however, the related spherical vector fields 

defined by ai = (^| — v'^)T -^UsU, a^+i = San induce exactly the mKdV hierarchy (and 
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this links our work with that of Chern-Tenenblat [34] on integrable fohations of surfaces 

with constant Gauss curvature). The two famiUes are related via cu^+i = Ti^io-n)- 

Yet another hierarchy of spherical evolution equations is obtained by simply re- 
stricting the even vector fields X2n to curves with v — vq. It is worth remarking that 
both the an and ujn can be expressed in terms of the restricted X2n- 



17] ^2n-2k 



k=0 
"■ / 1\k,,2k{ 1 



k\ 



M ^2n-2k 



k=0 

Here, {a)k = (a)(a + l)(a + 2) . . . (a + /c — 1) is the standard Pochhammer symbol. 

(v) Constant curvature preserving: The interest of the spherical evolution equa- 
tions discussed above is enhanced by the general connection between spherical curves and 
space curves, via the tangent indicatrix construction. To be brief, we mention only the 
following special case: a unit speed curve 7 on the unit sphere, with geodesic curvature 
u, is the derivative of a unit-speed space curve F, having constant curvature k = 1 and 
torsion r = u. It follows that we may regard the spherical equations of (iv) as evolution 
equations on space curves of constant curvature k = 1; moreover, in the case of the o"^ 
hierarchy, the torsion of F evolves according to mKdV. We note that it is also possible 
to translate the an into explicit (non-local) vector fields on the space curve level; e.g., ai 
corresponds to Ei = J ((^^2 — 1)-^ + '^sB ) ds, etc. 

Surface invariants. Given an evolving space curve, one can also consider the 
geometric invariants of the surface being swept out by the curve. For example, in case 
(iv), the curves are sweeping out a sphere, a surface of constant positive Gauss curvature. 
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We now discuss means of generating constant negative Gaussian curvature 
(CNGC) surfaces. These surfaces have long been a subject of interest in differential 
geometry, from the early work of Backlund [35], to the recent work of Sym [36] and 
Melko- Sterling [37], who have approached their study using some of the tools of modern 
soliton theory. Here, we indicate a complementary avenue of investigation: we discuss 
evolution equations on curves which sweep out CNGC surfaces. 

We begin by observing that the variation formulas (c),((i),(e) clearly simplify 
in the case that W is an eigenvector for TZ. Specifically, if 7 is a curve with constant 
torsion r = tq, one easily verifies that the vector field W = cos{9)T — sin{9)N satisfies 
71{W) = tqW, where = f^ K{u)du. Note that as a consequence of our variation 
formulas, the evolution equation jt — W induces an evolution on curvature Ht — — < 
Nj7l'^{W) >= — < NjTqW > = TQsin{9), so 9 satisfies the sine-Gordon equation 

9st = —Gsin{9), where G = — Tq. 

This well-known soliton equation has long been associated with CNGC surfaces. Let us 
substantiate the claim that the swept out surface M is indeed CNGC, with curvature G. 
Observe, via our variation formulas, that W{t) = 0. Thus constant torsion is preserved 
during the evolution. Next, the normal i/ to the surface is just the binormal B of our 
evolving curve; hence the second fundamental form of the surface, applied to the tangent 
T of our curve 7, is just < -du{T), T > = < -B,, T > = <tN,T >= 0. Thus T is an 
asymptotic direction oi the surface M. By the Beltrami-Enneper theorem [38], the Gauss 
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curvature of the surface M along 7 is the negative of the square of the torsion of 7; this 
is just G. This dynamical prescription for generating CNGC surfaces, although closely 
related to standard discussions of the sine-Gordon equation [11], [16], seems to be new (it 
was recently and independently derived by Segur and McLachlan [39] in their study of 
evolution equations associated with surface motion in R^). 

Let's consider this evolution equation 7t = VF when the initial curve is a Hasimoto 
filament [22], a curve 7 with curvature k — 2asech(as) and torsion r = tq. One easily 
checks (with appropriate choice of antiderivative = J^ nds ) that 6ss = a^sin{6). Thus 
W = cos{9)T — sin{6)N = — a~^((-^| — a'^)T + k'N). Consequently, modulo scaling and 
slippage along the curve, W = Qi, one of our previous planar- like vector fields. 

Along 7, we have Oi = X2 — tqXi. Both X2 and Xi are Killing fields (= in- 
finitesimal isometries) along 7 [29]; hence the curvature and torsion functions remain 
unchanged (modulo translation) as the curve evolves. The comments of the previous 
paragraph therefore apply for all time; and the result is that the Hasimoto filament, 
evolving via the equation 7^ = Oi, sweeps out a CNGC surface. We contrast this result 
with an observation of Melko-Sterling: under the evolution 7^ = Xi, the Hasimoto fil- 
ament sweeps out a surface M which is a formal boundary point of the set of CNGC 
surfaces, although M itself is not CNGC. As our simple example suggests, there is a 
connection between the hierarchies of vector fields described in this paper, and the ge- 
ometry of CNGC surfaces - the vector fields Qn play a central role. Details will appear 
in a subsequent publication. 
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